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I. INTRODUCTION 



A. THE MAIN IDEA 

Image processing by nonoptical means has received exten- 
sive attention in the last few years. Several books and many 
papers have been published that have established nonoptical 
image processing as a viable area of research. A large 
portion of this research emphasizes the linear processing 
of images for two main reasons: 1) Many image processing 

tasks are linear in nature. These tasks include image enhance- 
ment, image restoration, picture coding, linear pattern 
recognition, and TV bandwidth reduction. 2) There are many 
known linear techniques that may be brought to bear in the 
treatment of linear image processing. These techniques include 
transform theory, matrix theory, filtering, signal modeling, 
etc. Several common operations involved in image processing 
include transfer function concepts, partial difference (recur- 
sive) equations, and convolution summations. For example, 
Vander Lugt [Refs. 1,2] has presented an extensive development 
of linear optics based on transfer functions. The transfer 
functions relate the two-dimensional Fourier transform of an 
output image to that of the input image. Complex optical 
systems are easily described by combinations of transfer 
functions that correspond to individual components of the 
optical system. 
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Partial difference equations are used by Habibi [Ref. 3] 
to describe a model for estimating images corrupted by noise. 

The model corresponds to a two-dimensional extension of Kalman 
filters. Convolution summations are discussed by Fryer and 
Richmond [Ref. 4] in work that involves simplifying a two- 
dimensional filter to a single dimensional filter. 

The time-discrete state-space model offers great 
utility in the formulation and analysis of linear sys- 
tems. Linear systems that are described by transfer functions, 
difference equations, or convolution summations are formulated 
into a state-space representation. Once formulated, many 
known techniques may be applied to systematically analyze 
the model. Consequently, the state space model is a general 
and powerful tool that is used to unify the research and the 
study of time-discrete linear systems. 

This thesis develops the discrete model of Roesser [Ref. 

5] for linear image processing which closely parallels 
the well-known state space model for time-discrete systems. 
Because it is parallel, many of the concepts that are known 
for the temporal model may be carried over to the spatial 
model. This is done by generalizing from a single coordinate 
in time to two coordinates in space. The spatial model will 
hopefully have some of the same utility for the study of two- 
dimensional linear systems as the temporal model for one- 
dimensional linear systems [Ref. 3]. However, not all of the 
properties of one-dimensional systems carry over into the multi- 
dimensional case. 
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One of the fundamental problems involved with recursive 2- 
Dimensional systems is that the order of the system (recursive 
memory) is not the same as the number of initial conditions 
(boundary conditions) . In one-dimensional systems these are the 
same. Temporal systems are inherently nonanticipatory and are 
often treated as such for the sake of physical realizability 
in real time; whereas spatial systems do not have causality which 
is an inherent limitation. That is, an image processor may have 
right to left dependency as well as left to right dependency. 
Finally it is noted that stability criteria in one-dimensional 
recursive systems become much more difficult when carried over 
to the multidimensional case. 

Causality is built into the temporal state-space model if 
an initial state is assumed to be fully specified. In order to 
establish a close parallel for the spatial model, the same 
built-in causality will be intentionally assumed, despite the 
fact that causality is not necessary for physical realizability 
in real space. Such an image processor is said to be unilateral. 
If the constraint of causality is removed, then the image proces- 
sor is said to be bilateral [Ref. 5] . Concepts that are 
developed in this thesis for the latter case are: 

1) Formulation of the state space model of Roesser. [Ref. 5] 

2) The definition of state transition matrix. 

3) A resulting computer program based on the above model. 

4) An investigation of the class of 2 -Dimensional transfer 
functions defined by this model. 

5) Derivation of a general response formula. 
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6) Extension of Roesser's model of state variable equations 
to encompass a larger class of transfer functions. 

7) Adaptation of the 1-D "SSPACK" program to produce 2-D data. 

B. STATE SPACE REPRESENTATION 

Toward the end of the 1950s, the concept of representing 
a discrete system by a set of first-order difference equations 
became a standard tool of the research engineer. These tech- 
niques have since become generally known as state-space repre- 
sentations. Such representations have become increasingly 
important during the intervening years because^ they often allow 
one to carry out a meaningful system design entirely in the 
discrete-time domain (in comparison to popular Z-transform 
methods) . That this is important follows basically from these 
factors : 

1. The system may be nonlinear so that transformation 
methods are not directly applicable. 

2. Time-domain concepts often give one a better insight 
into the analysis and synthesis of the system (fre- 
quently with the aid of a digital computer) . 

3. Cases in which the initial conditions are non-zero may 
be handled straightforwardly. 

A state space representation of a system differs from the 

conventional representation. In a conventional representation 

only the relationships between the input and output signals need 

be known. On the other hand, the state-space representation 

gives a total description of both the internal as well as the 

external signals of a system. 
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C. STATE-VARIABLE REALIZATIONS- -THE CONCEPT OF STATE — 

In 1-D linear systems theory and control theory, the 
concept of a filter state has played an important role. 

Basically the filter state at any point in time contains 
all the information necessary to compute the remainder of the 
filter output signal, given the input signal. One dimensional 
single-input, single-output filter realizations based on a 
state variable model can be written in the form: 

x (k+1) = Ax(k) + Bu (k) (I. la) 

y (k) = Cx (k) + Du (k) (I. lb) 

This form relates the input u(k) and the output y (k) through 
a state vector x(k) . The state vector evolves in time accord- 
ing to equation (I. la). The matrices A, B,and C and 1 xl matrix 
D govern the exact form of the input-output relationship. 

(In general these matrices may vary with the index (k) and 
the input and output signals may be vectors as well.) Quite 
often the components of the state vector are taken to be the 
constants of the Z ^ delay operators in a flowgraph represen- 
tation of the 1-D filter. 

A classic problem in state-variable theory representation is to 
find the matrices A, B,Cand D which will realize a particular system 
function H(z) with a minimum number of state variables. A 
similar approach may be taken to develop a 2-D state-variable 
model . 
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A 2-D discrete system may be defined as a mathematical 
abstraction which utilizes three types of variables to repre- 
sent or model the dynamics of a discrete-time process. The 
three variables are called the input, the output, and the 
state variable. The input variables u(i,j), serve as external 
forces which influence the dynamics or motion of the system. 

The output variables y ( i , j ) are the characteristic variables 
which are directly observable (measurable) by the external 
observer. The state variables x(i,j) characterize the internal 
dynamics of the system and are to be selected according to the 
following rule. 

These variables are formulated in such a manner that, if 
one knows the values of the present state variables x(i,j) along 
with the values of the input variables u(i,j) then the output 
variables y(i,j> and the next state variables x(i,j) are com- 
pletely determined. Moreover, the number of state variables 
used in a state-space representation must be minimized. A 
state-space representation may be visualized in block diagram 
form, as shown below. 




u«(U 



INTERNA L 

STATES 
X,(U) , X^i/i ) . . 




Figure 1.1 
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In Figure 1 . 1 , m-inputs , p-outputs and n-state variables are repre- 
sented. However, we willbe mainly interested in those systems which 
have one input (m = 1) and one output (p = 1) . It is important 
to note that the input and output variables appear external to 
the system, while the state variables are generally internal. 

The different input variables will be represented by the 
input vector u(i,j) where. 



u(i, j) 



u ]_ ( i , j ) 

u 2 (i , j ) 



the output vector y(i,j) where. 



y(i/ j) 



Y 1 (if j) 

V i,j) 



and the state vector x(i,j) where. 



x(i, j) 



x 1 (i, j) 

x 2 (i, j) 

x n U,j) 



For a given process the state space representation is not unique. 
However all such representations have one characteristic in 
common for a given system, namely the number of elements n is 
referred to as the order of the system. 
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II. ROESSER'S STATE- SPACE MODEL 



A. THE FRAMEWORK 

An image is a generalization of a temporal signal, in that 
it is defined over two spatial dimensions instead of a single 
temporal dimension. Consequently, two space coordinates i and 
j take the place of time, t. Also, two-state sets are intro- 
duced to replace the single-state set. The following defini- 
tions are made by the model: 

i An integer-valued vertical coordinate; 

j An integer-valued horizontal coordinate; 

{r} A set of n-|_ real vectors which convey information 
horizontally; 

{ S } A set of real vectors which convey information 
vertically ; 

{u} A set of m real vectors that act as inputs; 

{y} A set of p real vectors that act as outputs. 

A specific image processor is then defined as 6-tuple 

<{R} , {S} , {u} , {y } , f ,g> , 

where f is the next state function: 

f: {{R},{S},{u} -*• {{R},(S}} 

and y is the output function 

g: {{R},{S},{u}} - {y} . 
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Now since f and g are to be linear functions, they may be 



represented by the following matrix equations : 



R ( i+1 , j ) 
S(i, j+1) 
y ( i r j ) 



A 1 R(i,j) + A 2 S(i,j) + B 1 u(i,j) 

A 3 R(i,j) + A 4 S(i,j) + B 2 u(i,j) (II. 1) 

C 1 R(i,j) + C 2 S(i,j) + Du ( i , j ) i,j _> 0 



A-^, A 2 , A^ , A^ , B^ , B 2 , C^, C 2 , D are matrices of appropriate 
dimensions. Boundary conditions R ( 0 , j ) and S(i,0) and also the 
input u(i,j) are externally specified. In the next section a 
computational rule is obtained that uniquely determines the 
states R(i,j) and S(i,j) and also the output y(i,j) (for i,j _> 0) 
from the boundary conditions (such as all zero). The equations 
produce a set of output vectors from the input vectors. 

This formulation is general so that any discrete linear 
image process may be so represented. Notation is condensed 
somewhat by introducing the following matrices and vectors: 



A = 




T(i, j) 



R(i, j ) 
S(i, j) 



T' (i, j) 



T' (i, j) 
AT ( i , j ) 



C = [c 1 c 2 ] 

R ( i+1 , j ) 

S (i , j+1) 

+ Bu (i , j ) 



y ( i , j ) = CT ( i , j ) + Du ( i , j ) 



B. GENERAL RESPONSE FORMULA 

A state- transition matrix A is defined as follows: 
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A = 




Then exponentiation A 1 ' ^ is defined as, 



A i' j = A 1 , 0 A i_1,j + A 0 , 1 A i,j_1 



(i / j ) > (0,0) 



A°'° = I ; A 1 ' ^ = A 1 ' -1 ’ = 0 for j >1, i > 1 



Examination of this definition bears out that it is an effective 
recursive definition of A 1 '- 1 for integer values of i and j such that 
either i > 0 or j > 0 or (i,j) = (0,0). It parallels the defini- 
tion of the time-discrete state-transition matrix A*" = A A fc ^ . 

It now remains to be shown that this state transition matrix 
A 1 '- 1 may be used in expressions for the response of the model 
in terms of the inputs and boundary conditions. The term 
boundary conditions is used here to refer to the states along 
the edges of the model. Specifically, the set of boundary condi- 
tions consist of R(0,j) for j >_ 0 and S(i,0) for i ^ 0. 



C. CHARACTERISTIC FUNCTION OF A MATRIX 

If the primary inputs and outputs are dropped in the model 
equations (II. 1), a representation arises for the state behavior 
of the system having the form 



R ( i+1 , j ) = A 1 R(i,j) + A 2 S(i,j) 
S ( i , j +1) = A^R ( i , j ) + A 4 S(i,j) 



(II. 2) 



These equations are useful in the development of a form for a 
two-dimensional characteristic matrix of A. Operators are 
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first introduced that advance a particular coordinate of 
their operand. 

Definition: Let E be an operator that has the effect of ad- 

vancing the vertical coordinate or the first subscipt of the 
function upon which it is operating. Likewise, let F be an opera- 
tor that has the effect of advancing the horizontal coordinate or 
second subscript of the function upon which it is operating. 

The effect of these operators on the state vectors is: 

R ( i+1 , j ) = ER ( i , j ) 

S ( i , j+1) = FS ( i , j ) 



The state equations can be rewritten using these advance 
operators . 



(EI-A 1 )R(i f j) - A 2 S(i,j) = 0 

-A 3 R(i,j) + (FI-A 4 )S(i,j) = 0 



These equations are equivalently represented in the overall 
matrix form. 



(EI-A 1 ) 




-A, 



(FI-A 4 ) 



T(i, j) 



0 



The above equation represents a system of homogeneous equations in 
the elements of T(i,j). If the system is to have a non-trivial 
solution for T(i,j) then the transformation represented by the 
matrix must be singular. The above matrix is said to be the two- 
dimensonal characteristic matrix of the partitioned matrix A, where 
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The characteristic matrix of A is denoted cm (A) and may be 
represented as 

cm (A) = El 1,0 + FI 0,1 - A 



where , 



i ! o 


c 


0 1 0 
I 


1 o 1 

,iT 


and I = 





Now since cm (A) must be singular, its determinant must be equal 
to zero. | cm (A) | = 0. If E and F are replaced in the above by 
general indeterminates x and y respectively, the result is an 
expression called the two-dimensional characteristic equation 
for A. The determinant of cm (A), and x and y replacing E and 
F, is called the two-dimensional characteristic function of the 
matrix and is denoted by 

| cm (A) | = f (x,y) = 0 

f(x,y) will be a monic polynomial in x and y with degree n^ 
in x, and degree n 2 in y, where n^ is the dimension of R and 
n 2 is the dimension of S. f(x,y) has the form 

f(x,y) = l l a. • x 1 y^ 

(0,0) < (i, j) < ( ni ,n 2 ) 1 ' J 



where a. • denotes elements of A and a „ = 1. 

i,D n l ,n 2 
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D. CIRCUIT ELEMENTS AND THEIR REALIZATION 

Let us consider the single 2-D HR filter transfer function 
given by: 



H ( z ]_ ' z 2 ) 



b 00 +b 10 z i 1+b 01 z 2 1+b ll z l l2 2 1+b 21 z i 2z 2 1 



-2 -1 



1 a 10 Z l -a 01 Z 2 a li z i Z 2 a 21 Z l Z 2 



B (z 1 , z ) 
1-A(z lf z ) 



A simple block diagram for follows. 






B (2., 20 



y(U) 



W (ijJ 



A U./ZO 



Figure 2.1 



The input signal u(i,j) flows through a filter corresponding 
to the numerator transfer function B(z^,z^). The resulting 
signal is added to the signal -w ( i , j ) to produce the output 
signal y(i,j). The denominator transfer function 1-A(z^,z^) 
is realized by the feedback loop containing A(z^,z^). 

Since we are dealing with two dimensions, there are two 
fundamental shift operators which may occur along a signal flow 
path, the horizontal shift operator indicated by z^ and the 
vertical shift indicated by z~^ [we shall omit from consideration 

i 

the inverse shift operators z-^ and z^\ . In most cases of prac- 
tical interest they can be eliminated by multiplying both the 
numerator and denominator polynomials of H(z^,z 2 ) by the appro- 
priate powers of z^- and z ^~ . Let us look at a signal flowgraph 
representing the numerator polynomial : 
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B(z 1 ,z 2 ) 



= b 



00 



+ b 10 2 l 1 + b 01 z 2 X + b ll z l lz 2 1 



^ ~2 -1 
+ b 21 Z l Z 2 



which is shown in Figure 2.2 below. 




Note the chain of two z^ b operators descending on the left and 

the single z^ b operator ascending on the right. The nodes 

along these two vertical paths are connected by branches with 

the appropriate gains. If we label the nodes in both z^ b 

chains and the chain 0,1,2 and so on, from the top down, 

the ith node in the z” 1 chain is connected to the jth node 

in the z_^- chain by a branch with a gain factor of b. .. 

^ i J 

Similarly the signal flowgraph for the polynomial A{z^,z^) 
is shown in Figure 2.3. 
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Since there is no term, there is no direct connection be- 

tween the input and output nodes of this signal flowgraph. Thus 
any path from the input node to the output node will encounter 
at least one z^ b or z^~ shift operator. 

At this point it is appropriate to discuss realizations for 
the two shift operators z^ b and z ^~ . At their simplest level, 
the shift operators merely select the "previous" S-tuple value 
in the horizontal or vertical direction. When the input to a 
Zj^ operator is the S-tuple u(i,j) the output will be R(i-l,j). 
Similarly for a z^~ operator the output will be S(i,j-1) when 
the input is R(i,j) or S(i,j). Consequently a realization of 
either shift operator must embody the appropriate amount of memory 
to retain the "previous" S-tuple in the appropriate direction. 

Interestingly enough, in the more general case where the numer- 
ator and denominator polynomials are considered jointly, the state 
variable realizations based on conventional signal flowgraphs may 
not be minimal in the sense that the transfer furnction can be 
realized with fewer coefficients. Consider, 



H ( z ]_ ' z 2 ) = 



b 10 Z l +b 01 Z 2 +b ll Z l z 2 
: =i =1 -i -i 

1-a 10 z l a 01 z 2 a ll 1 Z 2 



(II. 3) 



The corresponding signal flow representation is shown in Figure 
2.4 below: 




21 



E. ANALYSIS OF ROESSER'S MODEL 



Recalling from page 14 the equations of the model are: 

R ( i+1 , j ) = A^R ( i , j ) + A 2 S(i,j) + 

S ( i , j +1) = A 3 R(i,j) + A 4 S(i,j) + B 2 u(i,j) 
y(i,j) = ^Rfijj) + C 2 S(i,j) + Du ( i , j ) 

A^ , A 2 , A^ , A^ , B^, B 2 , C^, C 2 , D are scalars or matrices of 



appropriate dimensions. 



R ( i+1 / j ) 
S(i, j+1) 



A 1 A 2 




R(i, j) 


+ 


B i 


A 3 A 4 




S (i, j ) 




_ B 2_ 






— , 







(II. 4) 



y (if j) 



[ c i c 2 ] 



R(if j ) 

s (i , j ) 



+ Du (i , j ) 



(II. 5) 



R ( i+1 r j ) = A^R ( i , j ) + A 2 S(i,j) + B 1 u(i,j) 

S ( i+1 , j ) = A 3 R(i,j) + A 4 S(i,j) + B 2 u(i,j) 

And taking Z transforms: 

Z 1 R ( z i ' z 2 ) = AfR ( / z 2 ) + Z 2 S(z 1 ,z 2 ) + B 1 u(z 1 ,z 2 ) 
Z 2 S ( z ]_ ' z 2 ) = A 3 R(z 1 ,z 2 ) + A 4 S(z 1 ,z 2 ) + B 2 u(z 1 ,z 2 ) 



y(z 1 fZ 2 ) = [C ± c 2 ] 



R ( z i ' z 2 ) 

s (z ± , z ) 



+ Du ( z-j^ » z 2 ) 



(II. 6) 
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or 




z 1 R(z 1 ,Z 2 ) ” A 1 R(z r z 2 ) - a 2 S(z 1 ,z 2 .) 


= B 1 u(z 1 ,z 2 ) 


Z 2 S(Z 1 ,Z 2 } ' A 3 R (2 ] _/ z 2 ) ~ a 4 s ( z 1 ' z 2 ) 


= B 2 u(z v z 2 ) 


or 





R ( Z i ' Z 2 ) [z l -A 1 J ~ A 2 S(z 1 ,z 2 ) = B 1 u(z 1 ,z 2 ) 



R(z 1 ,z 2 ) [-a.,] - [z,-A]S(z., ,z 0 ) = B 0 u(z.,,z.,] 



1 ' 2 ‘ 



2 12' 



or 



or 



z r A i 

“-A “ 



-A, 



2 2 ~ a 



R(z 1 ,z 2 ) 
S (z^z ) 



B. 



B. 



u( z i,z 2 ) 



z il 0 
~o ■ ; ~z~ 

I 



A 1 . A 2 

- - J . 

A 3 ; a 4 

I 



^(Z^/Z^) 

sT z i , z 2 ) 









B 1 




B 2 







u(z 1 ,z 2 ) 



Where Z^ = z^I and Z 2 = z 2 I, and 



R(z 1 /Z 2 ) 
S (z 1 ,z 2 ) 



z‘, 0 

— 1 - 

0 ‘X 

I 2 



A i ; a 2 

A 3~ ! _A 4 



-1 



B. 



B. 



u (z., z 2 ) 



and after substitution in Equation (II. 6) 



y ( z 1 , z 2 ) = [c x c 2 ] 



-l 



| 








— 


z x 1 o 




J 

A 1 ' L A 2 




B 1 


~°T Z 2 




A 3 \ A 4 




B 2_ 



u (z- L , z 2 ) 



+ Du(z- l /Z 2 ) 
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or 



H ( z ]_ ' z 2 ) 



Y ( Z 1 ' Z 2 ) 

u(z x , z 2 ) 



tC l C 2 ] 



+ D 



r— — 




- 




— l 


r 

1 


Z ! ! 0 




A 1 ! A 2 




_L 


B 1 


_ - 1 




— — | — 








0 i z 9 
1 




^3 * ^4 
1 






B 2 











(II. 7) 



The submatrix Z-^ is simply z^ times an identity matrix of the 
appropriate size. Similarly is z 2 times an, identity matrix. 
The objective of the state variable realization procedure is 
to find the matrices A, B, C, and D which yields an F(z^,z 2 ) 
that equals or approximates a desired system function Hfz^z^) . 

In essence, the equations of Roesser represent an implementation 
for which a design algorithm must be found. One choice for 
the state variables is the output signals from the shift 
operators . 

Thus R(i,j) is a vector containing the output signals 
from the z^~ operators and S(i,j) contains the output signals 
from the z^~ operators. (Note that the output signal of a 
shift operator signal path is not necessarily the same as the 
nodal signal at the node to which the signal path points.) 

If a state variable corresponds to the output of a shift operator, 
the next value of that state variable must correspond to the 
input of the shift operator. To obtain the submatrices A^, 

A 2 , A^ , A^ in equations of Roesser, we write the input signal 
of each shift operator in terms of the outputs of all the 



2 4 



shift operators, taking care to include all shift-free paths 
from ouput to input (see the following flowgraph) . 




Figure 2.5 



Expanding the form of Equation II-7, page 24, yields: 



K(z l ' z 2 ) 



1 




l 

A x i A 
i_ i 


< — i 
i 


B ! 






— 3 * 







A 1 = j; adj A 

det A J 
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-1 



[Cl C 2 ] 



( 2 r A x ) ( 2 2*^4 ) "^ 2^3 



i 




“ 


Z _ -A . 1 A. 




B n 


f 

» 

1 

“I 

^1 

1 

N 

1 






A _ 1 Z. ~ A . 




B« 


3 ,11 




2 



IC-l C 2 i 



<VV B 1 



A 2 B 2 



( 2 2 -A- l ) (A 2 -Z 4 )-A 2 A 



a 3 b i 



( 2 1 -A 1 ) (Z 2 -A 4 ) -A 2 A 3 



(Z 1- A 1 )B 2 



(Zj^-Aj^) (Z 2 -A 4 ) -AjA 



( Z 1 ~A 1 ) (Z 2 -A 4 ) -A 2 A 



[Cl c 2 ] 



(Z 2' A 4 )B 1 + A 2 B 2 



(z 2 -Ai> u 2 -a 4 ) -a 2 a 3 
A 3 E 1 + (Z 1- A 1 )B 2 



‘W (Z 2 ' A 4 , - A 2 A 3 



or 



H (z lf z 2 ) 



C l( Z 2- A 4 )B l + C 1 A 2 B 2 + C 2 A 3 B 1 + C 2 (Z 1 _A 1 )B 2 

(Z 2 -A 4 )-A 2 A 3 



or 



h( 2 i ' z 2 } = 



C 1 B 1 Z 2~ C 1 B 1 A 4 +C 1 A 2 B 2 +C 2 A 3 B 1 +C 2 B 2 Z 1 C 2 B 2 A 1 
Z 1 Z 2 ~ A 4 Z 1 “ A 1 Z 2 ” A 2 A 3 +A 1 A 4 



( c 1 A 2 B 2+ c 2 A 3 B i“ C 2 B 2 A 1 -C i B i A 4^ + ^ C 2 B 2 Z 1 +C 1 B 1 Z 2 ^ 



( a 1 A 4 a 2 a 3 ) a 4 Z i “ A i Z 2 +Z 1 Z 2 



(II. 8) 



Equating equation (II. 8) with (II. 3) on page 21 yields 
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( C 1 A 2 B 2 +C 2 A 2 B 1~ C 2 B 2 A 1~ C 1 B 1 A 4 J +C 2 B 2 Z 1 +C 1 B 1 Z 2 
(A 1 A 4~ A 2 A 3 ) ~ A 4*l “ A 1^2 + Z 1 Z 2 



b 10 1 +b 01 Z 2 +b ll Z l z 2 

i ,"1 “1 _1 _1 

1 a 10 1 a 0l 2 a ll Z l Z 2 



For this example, Z^ = z^, Z 2 = 


z 2 ' 


all of the 


coefficients on 


the left hand side are scalars. 


Equation terms 


of equal powers 


of z^ and z 2> 








C 1 A 2 B 2 +C 2 A 3 B 1 -C 2 B 2 A 1 _C 1 B 1 A 4 


= 


b ll " 0 ' 


(II. 9) 


C 2 B 2 


= 


b 10 


(II . 10) 


C 1 B 1 


= 


b oi 


(II. 11) 


A 1 A 4 " A 2 A 3 


= 


1 


(11.12) 


A 4 


= 


a io 


(11.13) 


A 1 


= 


a oi 


(11.14) 


a ll 


= 


-1 


(11.15) 


From these equations, assuming 


that 


B 1 = B 2 = 1 


, it follows that: 


i — 1 

u 


= 


b io 




C 2 


= 


b oi 




i — 1 
< 


= 


a oi 




A 4 


= 


a io 
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From Equation (11-12): 



A 1 A 4 " A 2 A 3 “ 1 - _a n 



' A 2 A 3 ~ a ll A 1 A 4 



A 2 A 3 



a n + a io a oi 



Let A 2 and take on particular values p and q respectively, 



A 3 ^ A 2 P 



or 



pq = a u + a 10 a Ql 



(II . 16) 



From Equation (II-6) : 



C 1 A 2 B 2 + C 2 A 3 B 1 C 2 B 2 A 1 C 1 B 1 A 4 b ll 



or. 



b 01 P + b 10 q b 10 a 01 b 0l a i0 b ll 0 (11.17) 



Substituting Equation (11.16) into Equation (11.17): 



a n +a io a oi 

b 01 q + b 10 q “ b 10 a 0l " b 0i a i0 ” b ll “ 0 



or 
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0 . (11.18) 



b-> ^1*3 “ n a m i n+foi i ) ] 



'10 



io oi oi io u ir 



(b oi a n +b oi a io a oi ) 



The results are just the same as in [Ref. 8]. After the com- 
parison between Roesser's model and the 2-D HR filter, 
described by Equation (II. 3), we have: 



A, = a 



01 



p ; pq = a 11 = a 10 a Q1 



'3 q? b 10 q (b 10 a 01 +b 0i a i0 +b ll )q+(b 0l a il +b 0i a i0 a 01 ) 



A. = a 



10 



C n = b 



01 



(II 



C 0 = b 



10 



B, = 1 



Bo = 1 



D = 0 



The foregoing equations relate the coefficients of the 2-D 
transfer function to the terms of the system matrices of the 
Roesser model. Equation (II. 1). 

Kung et al . [Ref. 2] have shown that the following state 
variable equations, which use only two shift operators, will 
also realize H(z^,z^). For the foregoing example, 



= 0 



• 19) 
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u(i, j) 



U(u) 



R(i+1, j ) 




a 0l ! P 




R(i, j) 




1 




= 









+ 




S(i, j+1) 




j ! a io 




S ( i , j ) 




1 



Y(i, j) 



[b 



10 



b oi ] 



R(i, j) 
S (i, j) 



or 



H ( z l' Z 2^ ^ b 10 b 01^ 



z i a io P 


i — 1 

I 


l 


Z 2" a 01 


- 


1 

-V. 



We can construct a signal flowgraph with only two shift 
operators. It is an equivalent figure to that on page 25. 

Kung et al . [Ref. 2] have also shown that state-variable 
realizations of the form of the equations above may be 
generalized for any system function E{z^,z^) which satisfies 
the following three conditions: 



°<.o 
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0, must be 



1) The constant term in the numerator, b^Q = 
zero . 

2) The largest powers of z -^ , in the numerator and 
denominator polynomials, must be equal, and 

3) The largest powers of z^~ in the numerator and 
denominator polynomials must be equal. 

There is one potential difficulty with state variable reali- 
zations of this type. The nonlinear equations defining p and 
q may result in complex values for these constants. For 
example, when b^ Q = b Q1 = 1, b^ = 0, a 1Q = a Q1 = 2 and 

a ll = we get P = q* = 2 ±j. 
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III. THE PROGRAM OF ROESSER'S EQUATIONS 

WITH SCALAR COEFFICIENTS (FIRST ORDER) 



A. AN EXAMPLE 

For a 4 x 4 data field the S and R matrices are indexed 
as follows: 



j j 



1,1 1,2 1,3 1,4 




1,1 1,2 1,3 1,4 


2,1 2, 2 2,3 2,4 






2,1 2, 2 2,3 2,4 






i 


" 


3,1 3, 2 3,3 3,4 


< 


f 


3,1 3,2 3,3 3,4 


4,1 4,2 4,3 4,4 




4,1 4,2 4,3 4,4 



S matrix R matrix 

For 4x4 Matrices 



The Initial Conditions are given by the values 
R ( 1 , 1 ) , R ( 2 , 1 ) , R ( 3 , 1) , R( 4 , 1) 

S(l,l) , S ( 1 , 2 ) , S ( 1 , 3 ) , S ( 1 , 4 ) 

The 2-D state variable equations can be written as: 



R( i+1 , j ) 
S ( i , j + 1) 



y(i/j) = 



= A 1 R(i,j) + A 2 S(i,j) + B 1 u(i,j) 
= A^R ( i / j ) + A 4 S(i,j) + B 2 u(i,j) 




R(i, j) 

S (i, j) 



The input 2-D field is taken to be, 



u(i,j) = 1, for i = j = 1 

= 0, otherwise. 
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The output data field is indexed as: 




B. THE 2-D FOURIER TRANSFORM 

The 2-D discrete Fourier transform Y(m,n) pf the output 
y(i,j) can be written as. 



M-l N-l - j 2 tt 

Y (m,n) = l l y ( £ , k) e 

£=0 k=0 



£m . „ kn 

a 

e 



or for convenience. 



Y(m,n) 



M 

l 

1=1 



N 

l y(£,k) 

k=l 



j 2 77 



(£-1) (m-l) 
M 



— j 2 TT 



(k-1) (n-l) 
N 



Y (m, n) : 2-D D.F.T. {y (i, j) } 

MxN: The dimension of the given data y(£,k) and D.F.T. 

Y(m,n) also. 

y(£,k) : Given data (The output as described above) . 

To develop the D.F.T. for two-dimensional signals we 
consider a finite area sequence y(£,k) which is zero outside 
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it is of area 



the interval 0 < Z < M-l, 0 < k < N-l, i.e. , 



(M,N) and construct the periodic sequence: 



yU,k) = y[ ( (£) ) M ( (k) ) N ) ] 



The original sequence y(£,k) is recovered by extracting one 
period of y(£,k), i.e.. 



We then define the discrete Fourier transform of y(£,k) to 
correspond to the Fourier series coefficients of y(£,k). 
However, just as we did with one-dimensional sequences, we 
will maintain the duality between the time and frequency 
domains by interpreting the D.F.T. coefficients to also be a 
finite 2-D sequence. Thus with Y(m,n) denoting the D.F.T. 
of y(£,k), we can write 



y U ,k) 



y u,k) R M , N U,k) 



1, 0 < i < M-l, 0 < . k < N-l 




0, otherwise 



£=0 k=0 




or 



M-l N-l 

y(£,k) = mN I Z Y ( m , n ) e 

m=0 y=0 
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or 



Y (m, n) 



M N (m - 1) - j 2TT- ( - k - 1) (n -l } 

l l y ± • ( ^ / k) e N e N 

1=1 k=l 



As an example, consider the case for M = N = 5. 
Given 2-D Data Sequence 



1.1 1,2 1,3 1,4 1,5 

2.1 2,2 2,3 2,4 2,5 

y(£,k) = 3,1 3,2 3,3 3,4 3,5 

i 4,1 4,2 4,3 4,4 4,5 

5.1 5,2 5,3 5,4 5,5 



Matrix 

M=5 



5x5 

N=5 



£=1,2, 3, 4, 5 k=l , 2 , 3 , 4 , 5 



Then , 



Y(l,l) = 
m=l,n=l 



y (1,1) + y (1,2) + y(l, 3) 
+y(2,l) + y (2,2) + y(2,3) 
+y (3,1) + y (3,2) + y(3,3) 
+y (4,1) + y (4,2) + y(4,3) 
+y (5,1) + y (5,2) + y (5,3) 



+ y (1,4) + y (1,5) 
+ y (2,4) + y (2,5) 
+ y (3,4) + y (3,5) 
+ y (4,4) + y (4,5) 
+ y(5,4) + y(5,5) 
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